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Abstract 

We present the numerical investigation of the fermionic two-body decays of tau sleptons fi,2 
and t sneutrino in the Minimal Supersymmetric Standard Model with complex parameters. In 
the analysis we particularly take into account the cosmological bounds imposed by WMAP data. 
We plot the CP-phase dependences for each fermionic two-body channel of fi,2 and r sneutrino 
and speculate about the branching ratios and total (two-body) decay widths. We find that the 
phase dependences of the decay widths of the third family sleptons are quite significant which can 
provide viable probes of additional CP sources. We also draw attention to the polarization of the 
final-state tau in the fi,2 decays. 

Key Words: CP-phase, sleptons decays, WMAP-allowed band. 



1 



1 Introduction 

The experimental HEP frontier is soon reaching TeV energies and most of the physicists expect that 
just there theoretically proposed Higgs bosons and superpartners are waiting to be discovered. There 
are many reasons to be so optimistic. First of all, in spite of its remarkable successes, the Standard 
Model has to be extended into a more complete theory which should solve the hierarchy problem and 
stabilize the Higgs boson mass against radiative corrections. The most attractive extension to realize 
these objectives is supersymmetry (SUSY) pQ. Its minimal version (MSSM) requires a non-standard 
Higgs sector [2 which introduces additional sources of CP-violation[3l [4] beyond the 8c km phase [5]. 
The plethora of CP-phases also influences the decays and mixings of B mesons (as well as D and K 
mesons). The present experiments at BABAR, Tevatron and KEK and the one to start at the LHC 
will be able to measure various decay channels to determine if there are supersymmetric sources of CP 
violation. In particular, CP-asymmetry and decay rate of B — > X s j form a good testing ground for 
low-energy supersymmetry with CP violation [6] . The above-mentioned additional CP-phases explain 
the cosmological baryon asymmetry of the universe and the lightest SUSY particle could be an excellent 
candidate for cold dark matter in the universe [7J |S] • 

In the case of exact supersymmetry, all scalar particles would have to have same masses with their 
associated SM partners. Since none of the superpartners has been discovered, supersymmetry must 
be broken. But in order to preserve the hierarchy problem solved the supersymmetry must be broken 
softly. This leads to a reasonable mass splittings between known particles and their superpartners, i.e. 
to the superpartners masses around 1 TeV. 

The precision experiments by Wilkinson Microwave Anisotropy Probe (WMAP) [9] have put the 
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following constraint on the relic density of cold dark matter \j_ 

0.0945 < QcDAih 2 < 0.1287 (1) 

Recently, in the light of this cosmological constraint an extensive analysis of the neutralino relic 
density in the presence of SUSY-CP phases has been given by Belanger et al. [TU] . 

Analyses of the decays of third generation scalar quarks [12] and scalar leptons [13] with complex 
SUSY parameters have been performed by Bartl et al. In this study we present the numerical in- 
vestigation of the fermionic two-body decays of third family sleptons in MSSM with complex SUSY 
parameters taking into account the cosmological bound imposed by WMAP data. Actually, we had 
performed some studies in this direction for squarks |14[ 115) incorporating all the existing bounds on 
the SUSY parameter space by utilizing the study by Belanger et al. [10] before. These investigations 
showed us that the effects of Mi and its phase <pu(i) on the decay widths of squarks are quite signif- 
icant. Now we consider third generation sleptons. Namely, we study the effect of Mi and its phase 
<Pu(i) 011 the decay widths of fi^ and v T . 

In the numerical calculations, although the SUSY parameters fi, Mi, M%, and Af are in general 
complex, we assume that fj,, Mi and Af are real, but Mi and its phase <Pu(i) take values on the 
WMAP- allowed bands given in Ref. [TDJ. These bands also satisfy the EDM bounds [16]. The 
experimental upper limits on the EDMs of electron, neutron and the 299 Hg and the 205 Tl atoms 
may impose constraints on the size of the SUSY CP-phases [T7] [IS]. However, these constraints are 
highly model dependent. This means that it is possible to suppress the EDMs without requiring 
the various SUSY CP-phases be small. For example, in the MSSM assuming strong cancellations 
between different contributions [19], the phase of \x is restricted to \(p^\ < 7r/10, but there is no such 
restriction on the phases of Mi and Af. In addition, we evaluate the parameter Mi via the relation 

x In our calculation, we have used WMAP-allowed bands in the plane M\ — ip which are based on 1st year data. Now 
the WMAP 3rd year data is also available but the new WMAP + SDSS combined value for relic density of dark 
matter does not change the numerical results in Ref. jlOj . namely the WMAP-allowed bands. See "Note added" section 
of Ref. [10]. 
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M2 = (3/5)|Afi|(tan6>n/)~ 2 which can be derived by assuming gaugino mass unification purely in the 
electroweak sector of MSSM. It is very important to insert the WMAP- allowed band in the plane 
Mi — ip into the numerical calculations instead of taking one fixed Mi value for all y-phases, because, 
for example, on the allowed band for /i = 200 GeV, Mi starts from 140 GeV for ip — and increasing 
monotonically it becomes 165 GeV for ip = tt. In Ref.[10] two WMAP-allowed band plots are given, 
one for /i = 200 GeV and the other for fi — 350 GeV. For both plots the other parameters are fixed 
to be tan/3 = 10, m H + = 1 TeV, Af — 1.2 TeV and ip^—ipA f —^- We here choose the masses for ti j2 
slcptons as m.f 2 =1000 GeV and 771^=750 GeV. These 2 values lead to a sneutrino mass m^ T —TAh 
GeV for M L < Mg. 

2 Tau Sleptons and Tau Sneutrino Masses, Mixing and Decay 
Widths 

2.1 Masses and mixing in slepton sector 

The superpartners of the SM fermions with left and right helicity are the left and right sfermions. In 
the case of tau slepton (stau) the left and right states are in general mixed. Therefore, the sfermion 
mass terms of the Lagrangian are described in the basis (fL,f R ) as [20l [21] 
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(3) 
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where m T , e r , and 9\y are the mass, electric charge, weak isospin of the r-lepton and the weak 
mixing angle, respectively, tan/3 = V2/V1 with Vi being the vacuum expectation values of the Higgs 



fields Hf, i = 1,2. The soft SUSY-breaking parameters M^, and A T involved in Eqs. (3-5) can 
be evaluated for our numerical calculations using the following relations: 
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The f mass eigenstates f\ and T2 can be obtained from the weak states tl and tr via the f-mixing 
matrix 
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where 



and 



ip f = arg[M 2 L ] = arg[A r - M *(tan/J)- 2 ^] 
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(10) 



J\Ml R \* + (ml-Miy J\Ml R \ 2 + (ml-Ml L )* 

One can easily get the following stau mass eigenvalues by diagonalizing the mass matrix in Eq. (2): 



< 2 = 2 ( M " + M ™ * V(M| X - M\ R Y + i\Ml R V 



The i> r appears only in the left state. Its mass is given by 
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(11) 



m~ T = Mi + - cos(2^)m : 



(12) 



Note that in this work we neglect CP-violation effects related to flavor change. Besides that the 
scalar mass matrices and trilinear scalar coupling parameters are assumed to be flavor diagonal. 
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2.2 Fermionic decay widths of and v T 

The lepton-slepton-chargino and lepton-slepton-neutralino Lagrangians have been first given in Ref. 
1. Here we use them in notations of Ref. 13: 



C vI . ± = gu{l%P R + k%P L )xjd t + gd{i%P R + k u l3 P L )xT^ + h.c. 



(13) 



and 



C-Uxo = glWijPR + K 3 P L )X% + h.c. 



(14) 



where u (u) stands for (s)neutrinos and d (d) stands for charged (s)leptons. We also borrow the 
formulas for the partial decay widths of h (k = n and v T ) into lepton-neutralino (or chargino) from 
Ref. 13. The partial decay width for the decay r, — > x? + T (^r) is expressed as 
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(16) 



where A r = ±| is the helicity of the outgoing r, k(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz), H s = 
[m~. - (m^o + m r ) 2 ] 1/2 and H p = [raf. - (m^o - ttv) 2 ] 1 / 2 . 
The explicit forms of the couplings, ajj, b\ 3 and £Jj, are 



(17) 



where 



f r 

\ h RjJ 



(18) 
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and 

flj = -^(N 32 +t^9 w N n ) 

h T Lj = {h T Rj T =Y T N* 3 . (19) 

The partial decay width of f, into the chargino, f, -4- + f r , is obtained by the replacements 

a ij ^ij'' ^ij "~ * 0' TO x° ~ ^ m x~' ™ T ~~ ^ ^ anc ^ ^ T ~ ^ — 2 m (15) and Eq. (16) with the couplings 
^ also given in Eq. (17) and Eq. (18). 

The width for the r-sneutrino decay v r — > X°jV T is obtained by the replacements ajj —¥ aj, i>T. — > 0, 
mfj — > Ti£; T , m T —¥ and A r — > — i in Eq. (15) and Eq. (16), and that for the decay z> r — > Xj~r(A r ) 
by the replacements a^- — > 6^ — > fcj, — > TOj> T and m^o — > m x +- The coupling are now 

a] = -^(timewNji - N j2 ), k» = Y T U* 2 , l] = -V jX . (20) 

Here, N is the 4x4 neutralino mixing matrix, U and V are 2x2 chargino mixing matrices and 
Y T = m T / (V2mw cos (3) is the r Yukawa coupling. 

In this work we contend with tree-level amplitudes as we aim at determining the phase-sensitivities 
of the decay rates, mainly. 

3 Tau-Slepton and Tau-Sneutrino Decays 

Here we present the dependences of the fi t 2 and v T two-body decay widths on the (fiu(i) f° r A* = 200 
GeV and for fi = 350 GeV. We also choose the values for the masses (m^, m? 2 , m ^±: m x ±> TO x?) = 
(750 GeV, 1000 GeV, 180 GeV, 336 GeV, 150 GeV) for n = 200 GeV and (m fl , m f2 , m-±, m-±, m^o) 
= (750 GeV, 1000 GeV, 340 GeV, 680 GeV, 290 GeV) for fi = 350 GeV. The mass values of these 
r-sleptons lead to a sneutrino mass TO£ T =745 GeV. Note that although the neutralino and chargino 
masses vary with <pun\, these variations are not large. Therefore, as a final state particle (i.e., on 
mass-shell), we have chosen fixed (average) mass values for charginos and neutralinos. For both sets 
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of values by calculating the and Mg, values corresponding to and m^, we plot the decay 
widths for > Mg and < Mg, separately. We plot the (^[/(^-dependences of the v r partial 
decay widths only for < Mg. In the case of > Mg, the phase dependences do not change, 
but decay widths take larger values. In our figures, we display the slepton decay widths for the both 
helicity states of the outgoing r (t£ and tr). 
M L > Mg for n = 200 GeV: 

In Figure 1(a) we show the partial decay widths of the channels fi — )■ x{v Tl f\ — > x^r, n — > 
Xi T L,i?: ^2 — > X\V T , t 2 — > x^^r and f 2 — > Xi t lm as a function of <£[/(i) for = 200 GeV. In these plots 
some dependences on the <fiu(i) phase are shown. In order to see these dependences more pronouncedly 
we now plot two channels separately; namely fi — > Xi T R (the variations in the cross section are not 
large) and f 2 — > Xi t l (the variations are really large) in Figures 4(a)-(b). Here we consider the case 
Mi > Mg, where fi is mainly fij-like and f 2 is mainly 7x-like (VJ xl —VJ 22 Ki 0). In this case, the decay 
processes whose initial and final state helicities are the same, f 2 — > x^^t, t 2 — > Xi^n ^2 — > Xi t l 
and fi — > Xi r fl) have large widths, whereas those with opposite helicities, f 2 — > Xi r fi; ^i ~~ xT^, 
fi — > x^^r and fi — > Xi T i, have small ones. The reason for these large and small widths can be traced 
to the couplings ajj, bjj, £Jj and aj, fcj, Because of H s rzH p (since 2 ^>m T ) we can express the 
decay widths of n -> x</"r(A r ) as T(f 4 -> x"r(A r = l/2))cx|6[/ and T(f 4 -> x"t(A t = -l/2))a|af-| 2 . 
For example, r(fi — > Xi T i) (r(T 2 — >• Xi t -r)) is suppressed because it is proportional to the term, \a[ x \ 
(|fe 21 |), which includes small Yukawa coupling (Y T ). On the other hand, r(fi — > Xi t r) is proportional 
to the square of b\ x which depends on the combination TZ\ X B\ x + 1Z\ 2 B\ 2 contributing largely from 
its second term. Similarly, since H s ~H p , the decay widths of fj — > xj + v t can be expressed as 
r(Tj — > xj + v T )ctHl \l\ 3 \ 2 . The decay widths of fi — > xT^tj n — > X2 are also suppressed due to 
the very small Yukawa coupling (£[ 1 «F T 7^[ 2 *J7i 2 , £]" 2 «y T 7\L[ 2 *?7 22 ). 

Note that the decay width r(fi — > Xi t r) is 90-110 times larger than r(fi — > Xi T i) and r(f 2 — > 
XiTl) is 10-30 times larger than T(f 2 — > Xi T it)- From Figure 1(a) one can see that the branching ratios 
for f 2 are roughly B(f 2 -> x 2 v T ) ■ B(f 2 -> XiTl) : B(f 2 Xi"^) : S(f 2 ->• x?Tr) ~ 6 : 1 : 0.5 : 
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Figure 1: Partial decay widths V of the fi^ decays for tan/3 = 10, A T = 1.2 TeV, <p a1 =<pa t =0, 
m fl = 750 GeV, m f2 = 1000 GeV and M~ L > fj, = 200 GeV in (a) and (jt = 350 GeV in (b). 



0.03. 

Although the <Pu(X) dependence of T(i2 Xi t l) (r(fx — > Xi t r)) stems only from the parame- 
ters | -Nil | and \Ni2\ (|-^n|), the phase dependence is quite pronounced. Similarly, the fu{\) phase 
dependence of T2 — > XiiXs)^ stemmed only from the <Pu{i) dependence of \U\\\ {\U2\\) parameter 
is also quite pronounced. The decay width T(t\ — > Xi t r) takes its maximum (minimum) value at 
Vi7(i)~lf" (yc/(i)~f ) ( see Figure 4(a)). This <pu(i) value also corresponds to maximum (minimum) 
value of \b\i | 2 . In a similar way, the width T(t2 — > Xi t l) and its parameter | | 2 takes their maximum 
(minimum) value at <Pu(i)~ n (Vi/(i) = 0) (see Figure 4(b)). Hence, we can say that the phase <Pu{i) 
dependence of |ay| 2 and \bjj\ 2 (|^L| 2 ) reflects the phase <pu(i) dependence of channels fj — > Xj r R,L 
(n -> Xjvr) directly. 
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Figure 2: Partial decay widths V of the fi,2 decays for tan/3 = 10, A T = 1.2 TeV, f^=fA T ~0, 
m fl = 750 GeV, m f2 = 1000 GeV and M~ L < fj, = 200 GeV in (a) and (jt = 350 GeV in (b). 



Mi > Mg /or (j, = 350 GeV : 

We give the same partial decay widths in Figure 1(b) for [i — 350 GeV (See also Figures 4(c)-(d)). 
Here, too, f\ is mainly f^-like and fi is mainly f^-like because we still keep the case > Mf,. For 
(j, = 350 GeV the WMAP-allowed band [10] takes place in larger Mi values (~ 305 - 325 GeV) leading 
to larger chargino and neutralino masses. This leads to smaller H 2 S (since H% oc [mf. — ™|J) and, as 
a result, smaller widths for f i 2 decays compared with those for /1 = 200 GeV. 

As can be seen from Figure 4(d), the ifiu(i) dependence of the decay T2 — > Xi t l is prominent such 
that the value of decay width at <Pu(i) =7r is about 2 times larger than that at ipjj^=0. The decay 
widths T(t2 — > xItr), T(fi — > X\V T \ T(fi — > %2. v t) and T(t± — > Xi t l) are suppressed because of 
the same reasons mentioned above. The decay width of the process fi — > X2 v r is the largest one 
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among the f 2 channels and the branching ratios are B(f 2 — > x 2 v T ) : B(f 2 — > Xi t l) : -5(f 2 — > x x i/ r ) : 
S(f 2 — > x\t r ) w 2.4 : 0.8 : 0.1 : 0.02. 
M z < Mg /or /z = 200 GeV: 

We give 7^2 and v T decay widths as a function of <£>[/(i) m Figure 2(a) and Figure 3(a) respectively 
(for /j, = 200 GeV). In Figures 4(e)-(f) we plot two of them separately whose CP-phase dependences 
are not clearly seen in Figure 2(a). They, too, show the significant dependences on CP-violation phase. 
In this subsection we consider the case < Mg, where fi is mainly fx-like and f 2 is mainly f^-like 
(JZ\ 2 =1Z 21 ~0) . The decay width of the process f\ — > x 2 ~ v t is the largest one among the f lj2 channels; 
its decay width increases from 3.55 GeV to 3.8 GeV monotonically as ipum increases from to w. 
In this case (M^ < Mg), the width r(fi — > x 2 v r) is not suppressed because its l\ 2 term does not 
include Yukawa coupling (i\ 2 ~ ^11*^21)- The phase dependence of f 2 — > Xi T R can be seen clearly 
in Figure 4(f); r(f 2 — » XiTr) takes its minimum and maximum values at and at ^[/(i)^^? 

respectively, because the parameter |65i| 2 reaches its minimum and maximum at these (fiu(i) values. 

The branching ratios for fi decays are roughly B(f\ — > X^t) : S(fi — > Xi t l) : S(fi — > X\ v t) '■ 
B(h x"tr) « 3.8 : 0.7 : 0.3 : 0.02. 

In Figure 3(a) we give i> r decay widths as a function of f° r A* = 200 GeV. The phase 

dependence is more significant for the decay channels v T — > xJ T io ^ t — > Xi~tz, an d £ T - ► Xi^t- 
Analogously to the neutralino decays of fi j2 ; because of H s taH p (since mp T >m T ) we can express the 
decay widths of £ T -> x+r(A T ) as T{v T x°r(A T = l/2))oc|fc|| 2 and T{v T -> x°r(A T = -l/2))oc|^| 2 . 
To be more specific, r(£> T — > Xi~Tr) (r(iv — > X^ T i?)) is suppressed because it is proportional to the 
term \k"\ (\k 2 \) which includes small Yukawa coupling (Y T ). Since H s =H p for neutrinos, the decay 
width of i> r — > Xi + can be expressed as r(i> r — > x? + ^lOocfff \a\ | 2 . The dependences of 

T(D T -> x/ T i) (r(i> T — > Xi^t)) stems from the <£[/(i) dependences of ( | o-i | ) parameter and this 
parameter is quite phase-dependent. 

The branching ratios for v T decays are roughly B(D T — > x 2 t l) '■ B(y T — > Xi t l) : -B(£v — > Xi^r) : 
B{v T £+7*) : B(j> t -> x+tr) « 3 : 1.3 : 0.1 : 0.01. 
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Figure 3: Partial decay widths T of the v T decays for tan/3 = 10, A T — 1.2 TeV, ififj,—fA T —0, — 750 
GeV, m f2 = 1000 GeV and ms T = 745 GeV; fi = 200 GeV in (a) and fx = 350 GeV in (b). 

M L < M & for (j, = 350 GeV : 

We present the dependences of the fi,2 and v T partial decay widths on <pu(i) m Figure 2(b) and 
Figure 3(b) (for /i = 350 GeV). In this case H^(Hp) takes smaller value because of the reason mentioned 
in the previous subsection and this leads to smaller widths for f\$ and v T . In Figures 4(g)-(h) we 
again plot two decay channels separately whose phase dependences are not clearly seen in Figure 
2(b). The dependence of the phase (fu(i) hi decays are similar to those in the case < Mg (for 
H = 200 GeV). 

Note that r(f 2 -► XirB)~ 90 T{f 2 -+ xW) and T(n -> x?r L )« 30 T(n -+ jgr fl ). 
The width T(i> T -> x°i + v T ) decreases as the phase increases from to n, showing a significant 
dependence on the phase. The branching ratios are roughly B(y T ~ > x\ v t) '■ B{p T — > Xi t l) '■ 
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Figure 4: (a)-(h) (fu(i) dependences of certain fi,2 decays for \i — 200 GeV and for [i = 350 GeV. 
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B{v T -> xtr L ) ■ B(D T -> x+t r ) : B{v T -> x+r fl ) ~ 0.6 : 0.2 : 0.15 : 0.09 : 0.0002. 

4 Discussion and Summary 

In this paper, we have presented the numerical investigation of the fermionic two-body decays of third 
family sleptons in the minimal supersymmetric standard model with complex parameters taking into 
account the cosmological bounds imposed by WMAP data. For this purpose, we have calculated 
numerically the decay widths of tau sleptons f\ t i and r sneutrino, paying particular attention to 
their dependence on the CP phase full)- We have found that some decay channels like r 2 — > X<2 v ti 

T2 -> Xl v T, h XlT L , f 2 -> Xl^R, h -> X2 V r, T X Xl^T, T\ -> Xl T L , V T -> X2 T L, V T -> x\ T L 

and v T — > x\ v t show considerable dependences on (pu(i) phase. These decay modes will be observable 
at a future e + e~ collider and LHC. Therefore they provide viable probes of CP violation beyond the 
simple CKM framework; moreover, they carry important information about the mechanism that brakes 
Supersymmetry. 

Besides that, f decay is important since it is the sole process where one can get information of the 
sfermion mixing and the neutralino mixing from the polarization of the final-state fermion |22) . Note 
indeed that for fi — 200 GeV and > Afg, the decay width r(fx — > Xi t r) 1S 90-110 times larger 
than r(fi — > Xi t l) and r(f 2 — > Xi t l) is 10-30 times larger than r(f 2 — > Xi T R.) since f\ (f 2 ) is mainly 
fjj-like ( rx-like). For [i = 200 GeV and < make only the interchange f\ <-> f 2 everywhere in 
the above-mentioned preceding sentence. For /x = 350 GeV the pattern expressed above remains more 
or less the same. 

The phase dependence of the fermionic two-body decay widths of f< and i/j stems directly from the 
parameters (Nij, Uij, Vij) of the chargino and neutralino sectors. The cosmological bounds imposed 
by WMAP data on the Mi parameter and its phase fuel) pl a Y an important role in taking their shapes 
of the phase dependences of these processes. 

In this study, we use the framework of R-parity conserving supersymmetric scenarios wherein the 
lightest supersymmetric particle (LSP) is a viable candidate for Cold Dark Matter (CDM). Other 
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than its relic density (observed by WMAP) little is known about the structure of CDM. But the 
recent astrophysical observations of the fluxes of high energy cosmic rays give information about the 
properties of CDM. In particular, recent results from Fermi LAT [23] indicate an excess of the electron 
plus positron flux at energies above 100 GeV. This also confirms the earlier results from ATIC [2"4] . 
On the other hand, PAMELA experiment |25j reports a prominent upturn in the positron fraction 
from 10-100 GeV, in contrast to what is expected from high-energy cosmic rays interacting with the 
interstellar medium. Although standard astrophysical sources such as pulsars and microquasars may 
be able to account for these anomalies, the positron excess at PAMELA and the electron plus positron 
flux of Fermi LAT have caused a lot of excitement being interpreted as decay/annihilation of Dark 
Matter. These unexpected results from PAMELA, ATIC and Fermi LAT experiments imply a new 
constrain on LSP that is CDM candidate: LSP must have not only the correct relic density found by 
WMAP but also correct decay/annihilation rates into electron-positron pairs. In the framework of the 
MSSM, a detailed analysis of decay of CDM that includes the observed cosmic ray anomalies has been 
given in Ref. 26 . 
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